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In the present paper an alternative theoretical interpretation to the generally assumed thermal emission models of 

the observed X-ray spectrum of isolated pulsar RX J0420. 0-5022 is presented. It is well known that the distribution 

function of relativistic particles is one-dimensional at the pulsar surface. However, cyclotron instability causes an 

appearance of transverse momenta of relativistic electrons, which as a result, start to radiate in the synchrotron 

regime. On the basis of the Vlasov's kinetic equation we study the process of the quasi-linear diffusion (QLD) 

00 ' developed by means of the cyclotron instability. This mechanism provides generation of optical and X-ray emission 

on the light cylinder lengthscales. The analysis of the three archival XMM-Newton observations of RX J0420. 0-5022 

is performed. Considering a different approach of the synchrotron emission theory, the spectral energy distribution 

. . is obtained that is in a good agreement with the observational data. A fit to the X-ray spectrum is performed using 

pL^ ' both the present synchrotron emission model spectrum absorbed by cold interstellar matter and generally assumed 

^ , absorbed black-body model. 

C^, Key words: pulsars: individual RX J0420. 0-5022, stars: magnetic fields, radiation mechanisms: non-thermal 

o 

t^ INTRODUCTION 

a 

The soft X-ray source RXJ0420. 0-5022 (hereafter RXJ0420) belongs to the group of seven radio-quiet, 
I ■ isolated neutron stars discovered in the ROSAT all-sky survey data and often called the X-ray Dim Isolated 
K^ I Neutron Stars (XDINSs, see [lOj for review). In the optical domain, only a few XDINSs have counterparts 
\^ , certified by their proper motion measurements and likely candidates are proposed for remaining, based on 
^^ ' their coincidence with the X-ray positions [17]. The X-ray spectra of XDINSs are well represented by an 
ff^ ' absorbed black-body {kT ~ 45 — lOOeV) which is emitted by the hotter part of the neutron star surface |10) . 
CN : Whereas, the cooler part of the XDINSs surface is assumed to be emitting in the optical domain [20]. The 
_il ■ formation of a non-uniform distribution of the surface temperature is still under debate, which is more likely 
(^ artificial and needs to be examined by convincing theory. The optical emission is also explained in terms of 
fT^ ' hydrogen layers of finely tuned thickness superposed on a condensed matter surface reprocessing the surface 
radiation [llj . or of non-thermal emission from particles in the star magnetosphere [[IBJ. 

Alternatively to somewhat contrived thermal emission models the observed properties of the faintest object 
in X-rays among the XDINSs RXJ0420 can be explained in the framework of the plasma emission model. 
This model has already been applied to explain the observed spectra of two other XDINSs (RX J1856.5- 
;h ' 3754 and RBS1774, see [H [5] for more details). The plasma emission model is based on the well-developed 
^_' theory of pulsars [141115) and suggests successful interpretation of the observational data. According to these 
works, in the electron-positron plasma of a pulsar magnetosphere the low frequency cyclotron modes, on 
the quasi-linear evolution stage create conditions for generation of the high energy synchrotron emission. 
Generally speaking, at the pulsar surface relativistic particles efficiently loose their perpendicular momenta 
via synchrotron emission in very strong magnetic fields and therefore, they very rapidly transit to their 
ground Landau state (pitch angles are vanishing). However, the cyclotron instability causes appearance of 
transverse momenta of relativistic particles in the outer parts of the pulsar magnetosphere. Therefore, the 
resonant electrons start to radiate in the synchrotron regime. 

We suppose that the observed X-ray emission of RXJ0420 is generated by the synchrotron mechanism 
at the light cylinder length-scales. In general the synchrotron radiation spectrum is considered to be a 
power-law, which is not consistent with the observational data from RXJ0420. The standard theory of the 
synchrotron radiation [2l [7] typically provides the power-law spectral energy distribution. Contrary to the 
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standard scenario, by taking into account the specifications produced due to the present emission model we 
obtain different spectral distribution, which can be successfully fitted with the measured X-ray spectrum of 
RXJ0420. 

In the present paper, we describe the plasma emission model and derive the synchrotron radiation spec- 
trum based on our scenario (Sec. 2), present the results of spectral analysis of re-processed XMM-Newton 
archival data from RXJ0420 (Sec. 3), and make conclusions (Sec. 4). 

EMISSION MODEL 

Any well known theory of pulsar emission suggests that, the observed radiation is generated due to pro- 
cesses taking place in an electron-positron plasma. It is generally assumed that the pulsar magnetosphere is 
filled by dense relativistic electron-positron plasma with an anisotropic one-dimensional distribution function 
(see Fig. 1 from [Ij) and consists of the following components: a bulk of plasma with an average Lorentz- 
factor 7p ~ 10^, a tail on the distribution function with 74 ~ 10^, and the primary beam with 7^, ~ 10^. The 
distribution function is one-dimensional and anisotropic and plasma becomes unstable, which might cause 
excitation of the plasma eigen modes. While considering the eigen-modes of the electron-positron plasma 
having small inclination angles with respect to the magnetic field, one has three branches, two of which are 
mixed longitudinal-transversal waves (/ti,2)- The high frequency branch on the diagram uj{k^^) begins with 
the Langmuir frequency and for longitudinal waves {k± = 0), Iti reduces to the pure longitudinal Langmuir 
mode. The low frequency branch, lt2, is similar to the Alfven wave. The third t mode, is the pure transversal 
wave, the electric component of which E* is perpendicular to the plane of the wave vector, and the magnetic 
field, (k, Bq). The vector of the electric field, E^^^''*^ is located in the plane (k, Bq). When k± = 0, the 
i-mode merges with the It waves and the corresponding spectra is given by 

0Jt^kc{l-6), ^- "^P 



2 0,3 



^^bV 



where k is the modulus of the wave vector, c is the speed of light, cOp = ^JA^:npe^ /m is the plasma frequency, 
e and m are the electron's charge and the rest mass, respectively, Up is the plasma density, ujb = eB/mc is 
the cyclotron frequency and B is the magnetic field induction. 

The main mechanism of wave generation in plasmas of the pulsar magnetosphere is the cyclotron insta- 
bility. The cyclotron resonance condition can be written as [1 2'| : 

U!r 
UJ - k\\V\\ - ka^Ua: + ^- = 0, (2) 

7r- 

where Ux = cV^'jr/pujB is the drift velocity of the particles due to curvature of the field lines with the 
curvature radius p. During the wave generation process, one also has a simultaneous feedback of these waves 
on the resonant electrons [24J. This mechanism is described by the QLD, leading to a diffusion of particles 
as along as across the magnetic field lines. Therefore, resonant particles acquire transverse momenta (pitch 
angles) and, as a result, start to radiate through the synchrotron mechanism. 

The wave excitation leads to redistribution process of the resonant particles via the QLD. The kinetic 
equation for the distribution function of the resonant electrons can be written as [5] : 

df^ d Id 

dt dp\\ '■ II J p^ dp± 

^ ^ ^P±D^,±^f'{p)]. (3) 



P± dp± { ' dp 

where 






P\\ 

are the transversal and longitudinal components of the synchrotron radiation reaction force and q^ 
2e2a;|/3c2. 
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Here -D_l,_l is the perpendicular diffusion coefficient, which can be defined as follows [1] 

D±,i. = ^JlEkf, (5) 

where \Ek\'^ is a density of electric energy in the waves and its value can be estimated from the expression 
\Ek\'^ ~ ?Tic^nf,7f,c/2a;, where oj is the frequency of the cyclotron waves. 

The transversal QLD increases the pitch-angle, whereas force Fj_ resists this process, leading to a sta- 
tionary state [df/dt = 0). The pitch-angles acquired by resonant electrons during the process of the QLD 
satisfies ip = p_l/p\\ ^ 1. Thus, one can assume that d/dp± >> d/dp\\. In this case the solution of Eq.(3) 
gives the distribution function of the resonant particles by their perpendicular momenta [3] 



where 

^ B-il V 32e6p3 

And for the mean value of the pitch angle we find -00 ~ PLq/p\\ — 10^^. Synchrotron emission is generated 
as the result of appearance of pitch angles. 

SYNCHROTRON RADIATION SPECTRUM 

To explain the observed X-ray emission of RXJ0420, let us assume that the resonant particles are the pri- 
mary beam electrons with 7^ r^ 10^, giving the synchrotron radiation in the soft X-ray domain. According 
to our emission scenario, the synchrotron radiation is generated as the result of acquirement of pitch an- 
gles by resonant particles, during the QLD stage of the cyclotron instability. As was shown in |16) . the 
cyclotron resonance condition (see Eq. (2)) is fulfilled at the light cylinder lengthscales. Consequently, the 
observed X-ray emission comes from the region near the light cylinder, where the geometry of the field lines 
is determined by the curvature drift instability excited at the same lengthscales |^19]. The curvature drift 
instability effectively rectifies the magnetic field lines (the curvature radius tends to infinity). Therefore, in 
the synchrotron emission generation region the field lines are practically straight and parallel to each other 
and one can assume that electrons with ^ ^ ipQ efficiently emit in the observer's direction. 
The synchrotron emission flux of the set of electrons in this case is given by (see [4J) 



/(px) = Cexp (J -^dpj] = Ce ^"-"0 ) , (6) 

^ vrVV 3mVSr , 



//i|(P||)-B^o- 



F, oc ' " ' ^ " ' ^ 



K5/3{z)dz 



dp\\- (8) 



Here ij|(p||) is the distribution function of the resonant electrons by their parallel momenta, e^ ~ 5 • 
W~^'^ Bipj'^keV is a photon energy of maximum of the synchrotron spectrum of a single electron and ^"5/3(2) 
is the Macdonald function. 

In order to find the synchrotron emission spectrum, one needs to know the behaviour of the distribution 
function, fAp^^)- For solving this problem, we consider the equation governing the evolution of / (pj) [1] 

dt dp^^ Vm2c27rV2'^^<"'iiy' ' ^ ^ 

Considering the quasi-stationary case from Eq. (9) we can find the redistributed distribution function of the 
resonant particles by their parallel momenta 

P\\ \Ek\ 
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On the other hand, the cyclotron noise is described by the equation 

^ = 2T,\E,\\ (11) 

where 

Tc = ^f\\iPr), (12) 

is the growth rate of the instabihty and from the resonance condition (2) it fohows that ku ~ us/cS^r- 
Combining Eqs. (9) and (11) it is easy to find that 0]: 

\Ek\'^pl'^, (13) 

here n denotes the index of the initial distribution function of the resonant electrons {fu oc P|7")- From 
Eqs. (10) and (13) it follows that /i|(p||) oc pn^ . As the emitting particles in our case are the primary beam 

electrons, nothing can be told about their initial distribution. We only know the scenario of creation of the 
primary beam electrons [8^, '3] which are extracted from the pulsar's surface via the electric field induced 
due to star rotation. To our knowledge there is no convincing theory that would predict the initial form of 
the distribution function of the beam electrons, which must be drastically dependent on the neutron star's 
surface properties and temperature. 

The frequency of the original waves, excited during the cyclotron resonance can be estimated from Eq. 
(2) as follows 

„ ~ 27r— ~ W^'^Hz. (14) 

o-Jb 

As we can see the frequency of the cyclotron modes comes in the same domain as the measured optical 
emission of XDINSs with the certified optical counterparts. Their spectra mostly follow the the Rayleigh- 
Jeans tail Fy cc v"^. On the other hand, the spectral distribution of the cyclotron modes is given by expression 
(13) and combining this equation with Eq. (14) we find 

Fj^oc |^fcpazv2"-3. (15) 

From Eq. (15) it follows that when n = 5/2 the spectral distribution of the cyclotron modes is coincident 
with the Rayleigh- Jeans function. And for the initial distribution of the beam electrons, as for their final 

distribution we find fu oc p,, and /ii oc pJ . We have used this distribution for the beam electrons to define 

the theoretical X-ray spectrum of RX J1856. 5-3754, which well fitted the measured one [5]. Although, for 
RXJ0420 it is not yet confirmed the detection of the optical counterpart, based on the similarity of XDINSs 
we assume that the initial distribution function of the beam electrons must be the same for all XDINSs. 
Thus, in place of integral (8) as for two other XDINSs investigated in previous works [3] [5] we get: 

F, oc 60-3 . e-(^A")\ (16) 

To find the values of the parameters b and tm, one should perform a spectral analysis by fitting the model 
spectrum absorbed by cold interstellar matter with the observed X-ray spectrum of RXJ0420. 

SPECTRAL ANALYSIS 

To obtain the X-ray spectra of RXJ0420 with the highest statistical quality we used the EPIC-pn data 
collected from three XMM-Newton observations between December 2002 and January 2003. The archival 
XMM-Newton data were processed with the Science Analysis Software (SAS) version 11.0. The X-ray spectra 
were grouped in spectral bins containing at least 30 photons. Subsequent spectral analysis was performed 
with XSPEC V12.7.0. The three EPIC-pn spectra of RXJ0420 were fit simultaneously and spectral analysis 
was limited to energies between 0.15 and 1.0 keV. Along with the plasma emission model proposed in the 
present work we also examined the black-body model, to compare the fit results. 
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Table 1: Model-Fit spectral parameters of RXJ0420. 0-5022 for combined fits to the EPIC-pn spectra obtained from 
individual observations in the energy interval 0.15 — 1.0 keV. 

Model Nh em b kT,f E.^ge r^dge x'(dof) 
(10^°CTO-'^) (keV) (eV) (keV) 

plasma l-OOtoJg 0.10 ±0.04 1.27 ±0.31 1.51(96) 

plasma*edge l.OOtoJy 0.10 ±0.05 1.27 ±0.42 0.30lJ^;j;i O.GSlJ^jJ^ 1.10(94) 

bbody l-Oltalg 43.3 ±1.2 1.47(97) 

bbody*edge I-Oltpja 46.3 ±1.8 0.3lj:°;g^ O-SQlpJl 1-10(95) 

THE BLACK-BODY MODEL 

For the spectral analysis of EPIC-pn data for RXJ0420 first we used an absorbed black-body model. We 
allowed only the black-body normalization to vary between the spectra of the individual observations and 
fitted the temperature and amount of interstellar matter as common parameters. We found that the black- 
body provided a reasonable fit, but with N^ ~ lO^^cm"^, that is not perfect. Thus, we set a lower limit 
for this parameter, which did not change drastically the fit quality. The resulting x^ = 1-47 for 97 degrees 
of freedom and for the column density we got Nh = (1.01 ± 0.19) x W^^cm^^. The best fit black-body 
temperature kT^ = 43eV appears to be the lowest value derived from any of the known XDINSs (see 
parameters in Table 1). 

THE SYNCHROTRON EMISSION MODEL 

The plasma emission model proposed in the present paper was recently applied to explain the X-ray spectra 
of RX J1856. 5-3754 and RBS1774, revealing good fit quality in both cases Ul [5]. We performed fitting of the 
model spectrum (Eq.(16)) absorbed by cold interstellar matter with the EPIC-pn spectra of RXJ0420. The 
best fit results are b = 1.27 and em = O.lkeV, corresponding x^ = 1-51 for 96 degrees of freedom. The column 
density Nh, Cm and b were treated as free parameters common to all three spectra. Only normalization was 
allowed to vary freely for different spectra independently, as in previous case when applying the black-body 
model. The fit results are listed in Table 1. 

DISCUSSION 

The spectral analysis of the pulse phase- averaged X-ray spectra of RXJ0420 shows that the fits quality in 
case of both models (black-body and plasma emission models) is not very good. From inspection of Fig.l 
one might consider the residuals around O.SkeV as an absorption feature. This feature, described as a broad 
(it = 70eV) Gaussian absorption line was noted and discussed by [9j. We add the absorption line at ~ O.SkeV 
to both models considered in the present work and re-fit the data. The fit quality improves in both cases, 
only if using an absorption edge model. The best-fitting energy of the edge is Eedge ~ O.SkeV, and the optical 
depth Tedge ~ 0.7 (see Tab.l). The resulting value of x^ reduces to 1.1 for both models. 

The nature of the spectral features discovered in X-ray spectra of XDINSs is not fully clarified as yet. 
According to the thermal emission scenarios the most likely interpretation of the absorption lines is that 
they appear due to proton cyclotron resonance. The proton cyclotron absorption line at ~ O.SkeV implies 
the magnetic field of Bcyc = 6.6 x lO^^G, which differs from the value of the dipolar magnetic field inferred 
from the timing measurements Bi^ip = 1.0 x lO^^G |10| . 

We suppose that existence of the absorption feature in X-ray spectra of RXJ0420 might be caused by 
wave damping at photon energies ~ O.SkeV, which takes place near the light cylinder. During the farther 
motion in the pulsar magnetosphere, the X-ray emission might come in the cyclotron damping range. If 
we assume that damping happens on the left slope of the distribution function of primary beam electrons 
(see Fig. 1 from [1]), then the photon energy of damped waves will be eo = {h/2TT)2ujBMb'4''^ — O.SkeV [3]. 
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Fig. 1: Combined EPIC-pn spectra of RX J0420.0-5022. Top left: The absorbed black-body model fit. Top right: 
The black-body model fit including an absorption edge at ^ 0.3keV. Bottom left: The synchrotron emission model fit 
absorbed by cold interstellar matter. Bottom right: Model fit including absorption edge at O.SkeV. 



Taking into account the shape of the distribution function of beam electrons, we interpret the large residuals 
around ~ O.SkeV as an absorption edge. 

Despite the fact that fit quality considerably improves when includin the absorption edge at 0.3keV for 
both emission models, the physical reahty of this feature is still uncertain. It might be caused by calibration 
uncertainties. A feature of possible similar nature was detected in EPIC-pn spectra of the much brighter 
prototypical object RXJ1856. 4-3754 and classified as remaining calibration problem by [10]. Thus, we agree 
that more data are necessary to finally prove or disprove the existence of this feature. 

According to the fit results (see Tab.l) the spectral analysis of the measured EPIC-pn X-ray spectra of 
RXJ0420 does not seem to be enough to distinguish between the black-body and the plasma emission models. 
Differently from the thermal radiation models, which appear to be somehow artificial, our scenario is based on 
a self-consistent theory. According to works of [21j and [22j due to the cascade processes of a pair creation the 
pulsar's magnetosphere is filled by electron-positron plasma with the anisotropic one-dimensional distribution 
function. The beam particles undergo drifting perpendicularly to the magnetic field due to the curvature, 
of the field lines. Both of these factors (the one-dimensionality of the distribution function and the drift 
of particles) might cause the generation of eigen modes in the electron-positron plasma if the condition of 
cyclotron resonance is fulfilled. The generated waves interact with the resonant electrons via the QLD, which 
leads to the diffusion of particles as along as across the magnetic field lines, and inevitably causes creation 
of pitch angles by resonant particles. Therefore, the resonant electrons start to radiate in the synchrotron 
regime. 

The estimations show that for the beam electrons with the average Lorentz-factor 7i, r^ lO"^, the syn- 
chrotron radiation comes in the same domain as the measured X-ray spectrum of RXJ0420. Differently 
from the standard theory of the synchrotron emission [7], which only provides the power-law spectrum, our 
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approach gives the possibihty to obtain different spectral energy distributions. In the standard theory of the 
synchrotron emission, it is supposed that the observed radiation is cohected from a large spacial region in 
various parts of which, the magnetic field is oriented randomly. Thus, it is supposed that along the line of 
sight of an observer the magnetic field directions are chaotic. In our case the emission comes from the region 
of the pulsar magnetosphere where the magnetic field lines are practically straight and parallel to each other. 
And differently from standard approach, we take into account the mechanism of creation of the pitch angles, 
which inevitably restricts their possible values. 
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